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WKOI  AIRPLANE  TURBTJLENCE  AND  FLUI  MEASUREKENTS 
O'NEILL,  NEBRASKA,  AUGUST  21  - 28,  1953 


I.  INTRODUCTION 

The  primary  purpose  of  this  report  is  to  distribute  data  obtained 
over  the  plains  of  Nebraska  to  other  participants  of  the  Great  Plains 
Turbulence  Field  Program  at  O'Neill.  Little  discussion  or  analysis  of  the 
data  can  be  given  until  more  data,  such  as  the  wind  profiles  obtained  by 
other  groups,  are  reduced  and  made  available.  A complete  analysis  vd.ll 
be  contained  in  the  final  report  of  the  Air  Force  Program, 

The  field  program  in  which  the  marine  meteorology  group  at 
Woods  Hole  was  invited  to  participate  was  conceived  and  organized  by 
members  of  the  Geophysics  Research  Directorate  of  tiie  Air  Force  Cambridge 
Research  Center.  They  and  their  contract  organizations  interested  in  the 
study  of  atmospheric  turtulence  set  up  their  equipment  in  an  east-vrest 
line  across  a flat  field  near  O'Neill,  Nebraska,  The  site  was  chosen 
for  its  extreme  flatness  and  the  predominance  of  southerly  •winds 
during  the  summer  months.  Figure  1 is  a photograph  vihich  shows  the 
observing  line  and  the  flatness  of  the  surrounding  land.  The  grovind 
stations  were  opera'ted  from  August  1,  1953  to  September  8,  1953. 

Woods  Hole  personnel  operated  only  between  August  21  and  August  28. 

During  this  interval  $6  horizontal  runs  were  made  between  UO  ft.  and 
7000  ft,  above  the  gro'und. 
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U.  THE  AIRPLANE  TECHNIQUE  OF  MEASUBINO  SHEARING  STRESSES,  ROOT-MEAN-SQUARE 

TURBULENCE  VELOCITT  COMPONENTS  AND  HEAT  FLUXES 

Only  bxrief  mention  will  be  made  here  of  the  airplane  method  of  neasxiring 
stresses  and  Telocity  conq>onents  since  the  technique  was  described  in  'UHOI 
Ref.  No.  53-52.  (About  20  copies  of  this  report  to  the  Office  of  Naval  Research 
are  still  available.)  It  is  planned  to  submit  a revised  copy  of  the  repoift 
to  a meteorological  journal  in  the  near  future. 

The  method  may  be  summarized  as  follows  t All  of  the  many  qaallflcatlons, 
limitations,  and  difficulties  of  the  system  will  be  passed  over  in  this  state- 
ment. Tilhen  an  airplane  flies  throu^  a turbulent  atmosphere  the  angle  of 
attack  continually  changes  thereby  modifying  the  lift  of  the  airplane.  If 
instruments  are  mounted  in  the  airplane  to  measure  thie  vertical  accelerations 
of  the  airplane,  the  airspeed,  and  the  attitude  of  the  airplane,  then  the 
vertical  velocity  of  any  gust  may  be  computed  by  placing  the  observed  values 
in  the  following  equation. 

* £l.l  Antj  - w„  (1) 

/ t yp  Vt*" 

Here  w^^  is  the  vertical  con?)onent  of  the  gust  iji  cn^/sec,  the  observed 
vertical  acceleration  of  the  airplane,  /=’  the  density  of  the  air,  the  true 
air  speed,  the  deviation  of  the  true  airspeed  from  a line  of  regz^sslon 

upon  time,  the  change  in  attitude  of  the  airplane,  t^  the  time  interval 

of  the  measurement,  and  w^  the  initial  vertical  velocity  of  the  airplane.  The 
coefficients  in  the  equation  depend  on  the  mass  of  the  airplane,  its  wing  f.roa 
and  other  considerations.  Values  of  w’  are  obtained  from  the  series  of  w^'s 
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by  eliminating  any  mean  value  of  and  finding  deviations  from  a line  of 
regression  of  upon  tima.  Tiiis  step  is  necessary  to  eliminate  the  effect 
of  any  mean  acceleration  of  the  airplane,  but  it  also  eliminates  the  con- 
tribution of  large  eddies  of  the  same  size  as  the  length  of  the  observing 
run  (1.5  km)  to  the  shearing  stress  and  root-mean- square  velocity  calculations. 

The  horizontal  conponent  of  the  gust  velocity  is  obtained  from  the  air- 
speed and  the  line  of  regression  of  the  airspeed.  Once  the  and  u^^  values 
are  obtained  they  majr  be  combined  to  form  the  shearing  stress,-^  w’u' , and  the 
root-mean-sq\iare  velocities,  and  Vu^.  A thermocouple  was  mounted 

on  the  airplane  so  that  deviations  of  the  temperature  could  be  combined  with 
the  w' s to  form  cp ^ w'T' , the  heat  flux. 

Probable  errors  due  to  reading,  calibrating,  airplane  deformation, 
instrumental  inaccuracies,  and  symmetry  assumptions  are  small.  The  errors 
arising  from  these  causes  based  on  l50  observational  points  amount  to  only 
+0.07  dynes  cra“2  for  the  sheartng  stress,  ^1.7  cm  sec"^  fori/  w>^  and 
+0.8  cm  sec”^  for  ^ u'  These  values  are  not,  however,  the  total  errors  of 

the  determinations.  Phugoid  motion  (a  vertical  oscillation  of  the  airplane 
of  about  20  seconds  period)  can  add  up  to  +20  cm  sec”^  to  the  root-mean -square 
velocities  for  an  extreme  case.  In  all  cases  presented  here  the  oscillation 
was  damped  so  that  it  is  believed  that  this  error  does  not  exceed  +5  cm  sec“^. 

On  August  25,  1953,  the  PBT  flew  past  a tower  on  which  Dr.  Soiani  had 
his  sonic  anemometer  operating  so  a direct  comparison  of  the  two  methods  of 
determining  w'  2 vas  made.  The  values  obtained  from  the  sonic  anemometer 
mounted  at  12  m were  72  and  68  cm  sec”^  for  the  first  and  second  run 
respectively.  The  PBY  values  were  58  and  60  cm  sec“^.  The  PBI  value  is 
based  on  a line  of  regression  against  tine  as  described  while  the  sonic 
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anemometer  is  based  on  a mean  valve  constant  with  time.  Vfhen  the  PBT  data 
is  reduced  in  a cojnparable  way  ?6  and  77  cm  sec"^  are  obtained  for  V w**. 

The  agreement  between  the  values  obtained  by  such  widely  differing 
methods  of  measurements  proves  that  the  airplane  method  is  reliable  to  at 
least  10$  accuracy  and  probably  better.  A frequency  diagram  (Figure  2)  of 
the  vertical  velocities  obtained  by  the  two  methods  is  presented  to  shew 
that  the  distributions  are  nearly  the  same. 

III.  EQUIPMENT  MOUNTED  IN  PBT  AIRCRAFT 

The  accelerations  of  the  airplane  were  measured  by  a Statham  electric 
accelerametar  mounted  near  the  center  of  gravity  of  the  airplane,  *1110  out- 
put of  this  strain  gauge  bridge  was  recorded  by  one  channel  of  a nine  charmel 
CEC  oscillograph.  For  the  air  speed,  the  dynamic  and  static  pressures  of 
two  pitot  static  tubes,  one  moimted  on  each  wing,  were  connected  to  a single 
Statham  differential  pressTire  gauge  to  give  an  average  value  over  the  wing- 
spread  of  the  air  speed.  To  record  the  attitude  of  the  airplane  a potenti- 
ometer circuit  of  an  automatic  pilot  control  was  connected  to  the  oscillograph, 
A reading  of  the  attitude  was  made  available  to  the  pilot  to  “fnable  him  to 
keep  the  airplane  at  constant  attitude.  The  thermopile  was  biased  by  two 
mercury  cells  and  the  resulting  signal  recorded  on  the  oscillograph.  The 
airplane's  altimeter,  airspeed  indicator  and  flux  gate  compass  were  used  for 
supplementary  observations. 

An  airplane  psychrograph  was  mounted  on  the  side  of  the  nose  turret 
to  give  dry-  ar  . wet-b\ilb  teirperatxrres  of  the  air.  Rod  thermistors  are  used 
in  a bridge  circuit  as  the  sensing  elements.  The  output  of  the  bridge  is 
recorded  on  an  L and  N Speedoraax  self-balancing  potentiometer. 
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IV.  FLIGHT  PHOCEDURE 

A standard  pattezm  was  developed  that  would  meet  the  reqtilrenents 

of  (1)  safet7  with  regard  to  both  ground  personnel  and  other  sounding  aircraft, 
(2)  non-interference  with  grottnd  observations,  and  non-interference  of  ground 
equipment  with  air  observations,  and  ( 3)  an  extended  horizontal  upwind  course 
over  terrain  comparable  to  that  upwind  of  the  gro\md  stations,  A course  that 
met  these  requirements  was  an  upwind  course  of  2-3  miles  extent  with  an  end 
point  200  or  300  yds  east  and  north  of  the  eastern  end  of  the  observing  line, 

A left  turn  was  made  at  this  point  and  the  airplane  climbed  to  the  starting 
point  of  the  course  for  another  run.  During  the  horizontal  runs  all  the 
equipment  was  operating.  The  psychrograph  vas  operated  alope  during  the 
200  ft  per  minute  climbs  to  the  next  run  level. 

Two  exceptions  to  the  above  flight  pattern  were  made  on  August  2?,  when 
comparison  tests  were  run  with  the  sonic  anemometer  of  the  University  of 
Wisconsin,  For  this  test  two  passes  were  made  past  the  instniment  mounted  at 
I4O  ft  on  a tower.  Actually,  the  fli^t  was  made  higher  than  the  instrument 
although  it  was  intended  to  be  at  the  same  level.  The  passes  xrere  made  v^)wind 
as  usual  and  the  instmiments  were  operating  both  prior  to  airival  at  the  tower 
and  for  a minute  or  so  after  its  pa,ssage. 

No  winds  were  measured  by  drift  meter  since  so  many  wind  observations 
were  made  by  the  gro\md  stations, 

V,  PRESENTATION  OF  REDUCED  DATA 

Both  the  graphical  and  tabular  forms  are  used  in  the  presentation  of 
the  reduced  data.  Certain  quantities  that  are  related  have  been  plotted  on 
the  same  height  diagram  to  clarify  their  relationship  and  allow  the  reader  to 


-6- 


get  a complete  picture  of  the  variatiOTi  of  the  quantities  throughout  the  layer. 
The  data  is  also  tabtilated  for  easy  refeience  and  extraction  of  exact  values. 

a)  Psychrograph  Data 

Values  of  the  potential  dry-bulb  temper atvire  and  the  mixing  ratio  are 
collected  in  Tp.ble  l.  Not  every  value  of  the  tempera t\ire  and  mixing  ratio 
is  given  in  the  table.  Rather,  only  enough  significant  levels  arc-  tabulated 
to  show  the  variation  with  hei^t  and  through  the  inversions.  It  will  be 
noted  that  all  available  values  have  been  plotted  in  the  diagrams.  The 
heights  given  in  the  taWe  are  heights  in  meters  above  the  ground  surface, 

b)  Turbiilenee  and  Flux  Data 

From  the  series  of  w^,  u^  and  T^  obtained,  the  following  quantities 
have  been  determined  and  entered  in  Table  II.  First,  after  an  identifying 
run  number  and  the  height  of  the  run  above  the  ground,  is  the  con^juted  value 
of  the  shearing  stress.  It  is  found  from  the  defining  equation  ■ -yc  w*u’ 
and  the  series  of  wj^  and  u^  values.  In  the  next  two  col\jmns  are  the  values 
of  l/  w>  2 and  ^ u— , Next  is  the  correlation  coefficient  betjween  the  w*  s and 
the  u*s.  In  the  next  column  is  the  heat  flux  computed  from  cp  w*T’ , In 
the  next  column  is  the  root-me an- square  deviation  of  the  dry-btCLb  temperature. 
In  the  last  column  a reliability  indicator  is  entered.  This  key  is  based 
entirely  on  the  change  in  attitude.  If  the  attitude  changed  by  less  than 
+1°  the  run  is  labelled,  R,  for  reliable.  If  the  attitude  changed  more 
than  +1°,  a U is  entered  indicating  an  unreliable  run. 

The  data  contained  in  Tables  I and  II  have  been  plotted  against  height 
in  seven  diagrams,  Figiurec  3-9.  On  the  left  is  the  potential  dry-bulb 
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Psychograph  Data  O’Neill,  Nebraska 

I 

TABLE  I 


A 

Aug.  21,  1953 
1112  to  lUiO  CST 


Height 

6d 

w 

m 

gny'Kg 

195 

30I1.7 

9.9 

2k0 

305.U 

9.8 

320 

305.0 

9.2 

500 

30li.9 

9.0 

61i0 

303.6 

9.I1 

770 

303.9 

9.1 

960 

30I1.I 

8.9 

1110 

30I1.2 

8.5 

1260 

30I1.2 

8.I1 

D 

Aug.  25,  1953 
1050  to  1150  CST 


9d 

w 

m 

oc 

gra/kg 

35 

307. U 

65 

307.2 

100 

307.0 

130 

307.0 

-3 

165 

307.3 

T> 

0) 

210 

307.1 

> 

360 

307.5 

•H 

h90 

308.7 

0 

660 

312.0 

825 

311.9 

0 

•H 

990 

312.3 

B 

Aug.  22,  1953 
Hi50  to  1<U0  CST 


Hei^t 

9d 

V 

... 

On 

tmlyro 
art  **o 

iti 

w 

35 

30T7.2 

7.8 

65 

307.2 

7.8 

1 An 

.^ww 

307.3 

7.7 

165 

307.3 

7.U 

285 

307.8 

7.5 

U80 

308..  2 

7.2 

700 

307.6 

7.a 

900 

307.5 

7.2 

1020 

307.6 

7.3 

1110 

307.1i 

7.U 

1260 

307.  U 

7.2 

Hi95 

307.0 

7.3 

E 

Aug.  25,  1953 
l600  to  1630  CST 


Height 

w 

m 

°C 

gm/kg 

180 

312.6 

360 

312.6 

li90 

312.5 

0 

660 

312.5 

825 

312. .5 

1000 

‘312.8 

♦Temperature  below 
180  m off  scsuLe  of 
psycho graph 


C 

Aug.  2ii,  1953 
15U7  to  l61i9  CST 


Height 

9d 

w 

m 

°c 

gm/kg 

35 

312.1 

11.3 

70 

312.3 

10.5 

135 

311.6 

10.7 

165 

312.1 

11.0 

360 

311.2 

10.8 

525 

312.1 

10.8 

675 

312.0 

10.5 

8h5 

311.7 

10.7 

1155 

311.6 

10.3 

1310 

312.2 

10.2 

1565 

311.9 

10.2 

1805 

311. 

10.1 

1950 

311.5 

10*0 

2070 

311.6 

9.9 

2130 

312.2 

8.8 

2215 

311.5 

6.7 

F 


Aug.  27, 

1953 

1110 

to  1211  CST 

Height 

Od 

w 

m 

°C 

gm/kg 

35 

308.2 

13.3 

65 

308,0 

13.5 

100 

308.5 

13.5 

130 

308.1 

13.5 

165 

308.1 

I3.U 

330 

308.6 

13.3 

li90 

308.5 

13.2 

660 

308.5 

13.1 

fi?o 

« 

juy.  V 

12.9 

990 

309.0 

12.6 

nUo 

309.0 

12.0 

1170 

311.  2 

10.3 

1200 

309.5 

11.8 

1260 

311.6 

10.0 

1365 

311.5 

9.6 

15I4O 

312.9 

8.5 

1670 

312.1 

8.3 
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Psychograph  Data  C 'Neill,  Nebraska 

TABIZ  I 
(continued) 


Q 


Aug.  28, 

1953 

10U3 

to  illi5  CST 

Height 

ed 

w 

m 

oc 

gm/kg 

3^ 

308.3 

17.1 

100 

307.8 

l6.0 

130 

307.6 

15.9 

180 

307.8 

l6.0 

330 

307.8 

15.7 

530 

308.2 

15.1 

660 

308.0 

llu7 

850 

308.0 

Ih.l 

990 

308.5 

13.7 

llliO 

312.  ii 

12.0 

1300 

312.  U 

11.5 

1500 

313.  U 

11.9 

l660 

3U.1 

10.6 
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teinperatupe.  All  values  measured  are  plotted  here  instead  of  the  few  signifi- 
cant ones  in  Table  I.  The  range  in  temperature  observed  on  any  hortzontal 
run  is  indicated  by  a line  with  an  "x"  at  either  end.  If  there  was  an  appreci- 
able change  in  altitude  the  line  is  drawn  sloping  to  indicate  the  altitude 
change.  Similarly,  the  mixing  ratios  are  plotted  both  as  points  and  range 
lines.  The  shearing  stress  and  the  heat  flux  have  been  plotted  with  a 
common  zero  line.  Different  symbols  and  scales  reduce  confusion  of  elements. 

On  the  ri^t,  and  /u'2  are  also  plotted  with  a common  zero  line  and 

different  symbols.  The  only  other  entry  on  the  diagram  is  a horizontal  line 
drawn  to  delineate  the  inversion  at  the  top  of  the  ground  layer.  This 
inversion  was  not  reached  in  all  soundings. 

Only  very  limited  discussion  can  be  made  of  the  data  \mtil  wind  profiles 
are  available.  Many  of  the  hei^t  variations  of  the  stress  values  seem  most 
peculiar  and  can  only  be  inteipreted  later  in  conjunction  with  the  measured 
winds.  It  is  hoped  that  such  huge  swings  as  the  stress  took  on  the  22nd, 

27th,  and  28th  will  be  e^lained  by  the  winds.  It  is  of  interest  to  note  that 
on  only  one  day  did  the  stress  have  a large,  positive  value  near  the  ground 
and  decrease  more  or  less  linearly  to  zero  near  the  ground  layer  inversion. 
This  distribution  of  stress  values  corresponds  to  the  variation  e:q}ected  if 
the  air  under  the  inversion  reacted  as  thou^  it  were  a liquid  confined  in  an 
open  channel,  an  assumption  frequently  made  in  dealing  with  the  atmosphere. 
This  observation  that  c^en  channel-like  flow  occurred  when  the  inversion  was 
relativeljr  low  may  give  a clue  as  to  the  role  of  the  inversion  height  in 
defining  a sin^e  flow  pattern  or  stratification  of  the  air  into  several 
distinct  flow  regions. 
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Values  of  - ^ w'u ' taken  in  or  near  an  inversion  should  not  be  interpreted 
as  a stress  or  fl\ix  of  momentiun.  This  is  true  since  the  airplane  frequentlj 
passes  from  one  air  layer  to  the  other.  Thus  correlation  between  w'  and  u* 
may  be  corap\ited  that  have  little  or  no  relation  to  the  actual  shear.  Similarly, 
heat  fluxes  ccinputsd  may  be  dependent  more  on  the  mtniher  of  times  the  airplane 
crosses  the  boundary  than  on  the  flow  of  heat. 


The  variation  of  the  values  of  'y  w'2  and  V u'2  may  be  commented  tpon 
briefly.  One  point  of  interest  is  that  the  'j/u*  ^ is  invariably  greater  than 
the  -yw'^  ac  the  30  m level,  but  at  either  the  6o  or  100  m level  thel^  w'2 
becomes  the  larger.  Above  the  100  m level  the  values  of  both  -^/w'^  and 
are  remarkably  constant  up  to  the  vicinity  of  the  inversion  where  they  decrease 
rapidly  to  a small  value.  When  values  of  the  coefficient  of  turbtilent  mass 
exchange  are  computed  it  may  be  possible  to  determine  the  relation  of  the 


V w*2  to  the  coefficient  without  relying  on  an  unobservable  mixing  length. 

No  spectrum  analyses  have  been  run  upon  the  present  series  of  data. 

The  limitations  of  the  system  are  such  that  only  a rather  narrow  spectral 
band  could  be  investigated.  The  smallest  averages  of  significance  may  be 
defined  as  either  10  or  20  m eddies  >diile  the  largest  are  of  the  order  1000  m. 
Spectral  analysis  within  these  limits  may  be  carried  out  with  assurance  of 
significance. 
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Legends  for  Figures 


Figure  1. 


Figure  2. 


Figure  3. 


Figure  li. 


Figure  5. 


Figure  6, 


Figure  7. 


Tbe  flatness  of  the  Nebraska  plains  can  be  Judged  from  this 
photograph  of  the  observing  line  at  O'Neill. 

The  frequency  distribution  of  the  turbulent  vertical  velocities 
as  measured  by  the  sonic  anemometer  and  the  airplane  are  com- 
pared. The  distributions  are  very  similar. 

Observations  of  stress,  root-inean-square  velocities,  heat 
flux,  potential  teit5>erature,  and  mixing  ratio  obtained  over 
O'Neill,  Aug.  21,  1953,  1012  to  lliiO  CST. 

Observations  of  stress,  root-mean-square  velocities,  heat 
flux,  potential  terape ratine,  and  mixing  ratio  obtained  over 
O'Neill,  Aug.  22,  1953,  lU50  to  l5U0  CST. 

Observations  of  stress,  root-mean-square  velocities,  heat 
flux,  potential  temperature,  and  mixing  ratio  obtained  over 
O'Neill,  Augc  2U,  1953,  15U7  to  l6U9  CST. 

Observations  of  stress,  root -mean-square  velocities,  heat 
flux,  potential  temperature,  and  mixing  ratio  obtained  over 
O'Neill,  Aug.  25,  1953,  iOU5  to  llii5  CST. 

Observations  of  stress,  root-mean-square  velocities,  heat 
flux,  potential  tenp^erature , and  mixing  ratio  obtained  over 
O'Neill,  Aug.  25,  1953,  l600  to  l650  CST. 


Legonds  for  Figures 
“Continued- 

Figure  8,  ObserrationB  of  stress,  root-mean-square  velocities,  heat  flux, 
potential  temperature,  and  mixing  ratio  obtained  over  O'Neill, 
Aug.  27,  1953,  1110  to  lUiO  CST. 

Figure  9.  Observations  of  stress,  root-mean-square  velocities,  heat 

flux,  potential  temperature,  and  mixing  ratio  obtained  over 
O'Neill,  Aug.  28,  1953,  10U3  to  1136  CST. 
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